Detecting stochastic multiresonance in neural networks via statistical complexity measure

This paper employs statistical complexity measure (SCM) to investigate the occurrence of stochastic multiresonance (SMR) induced by noise and time delay in small-world neural networks coupled with FitzHugh-Nagumo (FHN) neurons. Our findings reveal that SCM exhibits four local maxima at four optimal noise levels, providing evidence for the occurrence of quadruple stochastic resonances. When time delay \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau$$\end{document}τ is taken into account in the information transmission, under moderate noise levels, SCM shows several local maxima when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau = nT_{e}$$\end{document}τ=nTe with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document}n being a positive integer and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{e}$$\end{document}Te being the period of subthreshold signal. This indicates the appearance of delay-induced SMR at the multiples of the period of subthreshold signal. Intriguingly, at low noise levels, a strong coherence between time delay and neuronal firing dynamics emerges at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau = nT_{e} - 2$$\end{document}τ=nTe-2, as confirmed by a series of SCM maxima at these time delays. Furthermore, the study demonstrates that by adjusting the degrees and sizes of small-world networks, as well as the coupling strength, it is possible to optimize the strength of delay-induced SMR, thus maximizing the detection capability of subthreshold signal. The research results may provide us with an effective approach for understanding the role of time delay in signal detection and information transmission.


Detecting stochastic multiresonance in neural networks via statistical complexity measure
Yazhen Wu 1,2 & Zhongkui Sun 1* This paper employs statistical complexity measure (SCM) to investigate the occurrence of stochastic multiresonance (SMR) induced by noise and time delay in small-world neural networks coupled with FitzHugh-Nagumo (FHN) neurons.Our findings reveal that SCM exhibits four local maxima at four optimal noise levels, providing evidence for the occurrence of quadruple stochastic resonances.When time delay τ is taken into account in the information transmission, under moderate noise levels, SCM shows several local maxima when τ = nT e with n being a positive integer and T e being the period of subthreshold signal.This indicates the appearance of delay-induced SMR at the multiples of the period of subthreshold signal.Intriguingly, at low noise levels, a strong coherence between time delay and neuronal firing dynamics emerges at τ = nT e − 2 , as confirmed by a series of SCM maxima at these time delays.Furthermore, the study demonstrates that by adjusting the degrees and sizes of smallworld networks, as well as the coupling strength, it is possible to optimize the strength of delayinduced SMR, thus maximizing the detection capability of subthreshold signal.The research results may provide us with an effective approach for understanding the role of time delay in signal detection and information transmission.
It's well known that noise can lead to astonishingly constructive outcomes in nonlinear dynamical systems.Especially stochastic resonance (SR) 1 , as the most representative example of this fact, is a phenomenon that the response of a nonlinear system to a weak input periodic signal can be amplified by the assistance of optimal noise level.SR has been observed in various types of nonlinear systems, such as physical, chemical, mechanical, biological systems [2][3][4][5] .In the context of neural systems, SR has been documented in the caudal photoreceptors of crayfish and cricket 6,7 , predatory behavior of paddlefish 8 and mammalian brain 9,10 .Numerous experiments and theoretical studies have demonstrated that injecting appropriate noise into sensory neurons can significantly enhance their ability to process weak input signals 11,12 .Consequently, SR is considered as an important tool for detecting weak signals and gaining insights into information transmission within neural networks.
In recent decades, SR has been thoroughly explored in numerous neural network systems.In a comprehensive study by Perc et al., the effects of scale-free network topology on pacemaker-driven SR were thoroughly examined, revealing its remarkable efficacy in facilitating the transmission of weak localized rhythmic activity across the entire network 13 .Kwon et al. revealed that spatially correlated noise can significantly enhance the coherence resonance (CR) and reinforce the clustering structure of the network 14 .Furthermore, Liu et al. shed light on the potential of appropriately applied noise and electromagnetic induction in promoting SR and modulating the response of envelope modulation signals within neural networks 15 .Notably, a detail numerical investigation of SMR in a FHN neural network model was reported by Sun et al., who demonstrated that subthreshold signal could be detected successfully at more than one optimal noise level 16 .Together, these findings illustrate the substantial enhancement of information transmission or signal detection capabilities in neural networks through the application of noise-induced SR.
In fact, time delay is unavoidable during the process of information transmission in neural systems because of the limited speed at which action potentials propagate across axons.And it has been discovered that time delay has great influences on dynamical behaviors of coupled neural networks, such as enhancing or decreasing phase synchronization 17 , improving spatial-temporal order 18 , as well as inducing synchronization transition [19][20][21] .In particular, Franović et al. revealed that when there are interaction delays among neurons and they are influenced by external noise interactions, as well as with the excitability of individual neuron, synchronization clusters spontaneously emerge in the system 22 .Furthermore, time delay has been proved to play a crucial role in inducing SR in neural networks.For instance, Wang et al. discovered that time delay can induce the occurrence of SMR at every multiple of the pacemaker period in scale-free neural networks, facilitated by the presence of moderate noise 23 .Similarly, Liu et al. showed that when time delay is appropriately tuned in feed-forward-loop neuronal network motifs, multiple-SRs can appear 24 .Sun et al. reported that, regardless of whether the noise intensity was optimal or not, the detection capability of signals can be enhanced through the introduction of time delay in excitable neural networks 25 .Yu et al. revealed that the interaction between hybrid electrochemical synapses and time delay enhances the efficiency of signal transmission and promotes the detection capability of subthreshold signals, thereby strengthening the phenomenon of SR 26 .Recently, Tuo et al. demonstrated that the strength of delay-induced SMR can be controlled by adjusting the learning rate of synaptic plasticity in modular neural networks 27 .Yu et al. found that time delay can cause period variations in the signal transmission, resulting in SMR occurring at integer multiples of half the signal period in excitable neural networks 28 .Furthermore, they revealed that in delay-coupled neural networks, weak coupling can enhance the performance of signal responses within a broader range of delay windows, thus promoting the emergence of SMR induced by time delay 29 .Therefore, the investigation of the influence of time delay on SR in neural networks continues to be a prominent and popular research topic.
It is widely acknowledged that there are several traditional indicators, such as the signal-to-noise ratio (SNR) 30 , spectral power amplification 31 , Fourier coefficients 32 , and others 33 , which are employed to measure SR in neural networks.These indicators assist in the detection of weak signals and enhance our understanding of information transmission.However, there are limited studies focusing on SR in neural networks from the perspective of information theory, using SCM as a metric.Indeed, SCM is also an important quantifier to characterize SR, which has been profoundly investigated in low dimensional systems [34][35][36] .It has been illustrated that SCM could describe subtle features of internal systems.In a notable study, Sun et al. successfully showcased the efficacy of SCM in detecting double-SRs in fractional-order bistable systems, encompassing both strong and weak resonances 37 .In contrast, SNR can only capture the presence of strong resonances.This highlights that SCM is more effective and accurate than SNR in detecting weak signals within low-dimensional bistable systems.Building upon this, Wu et al. further expanded the application of SCM to measure SR behaviors in neural networks, validating its effectiveness 38 .However, their investigation did not consider the influence of time delay on SR in neural networks.
Inspired by the findings mentioned above, in order to reveal more subtle firing dynamical characteristics and delve deeper into the impact of time delay on weak signals detection in neural networks, it becomes crucial to investigate delay-induced SR using SCM as a tool in this study.Accordingly, the remainder of this paper is organized as follows: "Model and methodology" presents the mathematical model and introduces the utilization of SCM.Main results are exhibited in "Results".Finally, some conclusions are drawn in "Conclusion".

Mathematical model
The considered neural network is the ensemble of FHN neurons 39,40 , governed by the following equations Here x i and y i are separately the fast action potential variable and slow recovery variable of neuron i (i = 1, . . ., N) , with N being the total number of neurons in the network.ε represents the small-time scale ratio between the fast and slow variables.By fixing ε at 0.01 in this study, we are able to accurately describe the rapid evolution of membrane potential in the FHN neural network model and capture the dynamic characteristics of the system.g is the coupling strength.In this work, we consider a Watts-Strogatz (WS) small-world network with N neurons.To obtain it, a regular ring-like network is firstly constructed, where each unit is connected to its k nearest neighbors, then rewiring each edge randomly with the probability p .Based on the employed network, if neuron i is connected to neuron j , then J ij = J ji = 1 , otherwise J ij = J ji = J ii = 0 .τ is the time delay.The behavior of each individual neuron is determined by the system parameter a , where |a| < 1 signifies oscillatory behavior and |a| > 1 indicates excitable behavior.In this study, we assume that all units exhibit excitable behavior by fixing a value of a = 1.1 for every unit.A sin 2π T e t represents a subthreshold signal with an amplitude of A and a period of T e .That is, this signal cannot initiate neuronal firing without the presence of accompanying noise.Based on the research findings provided by Li et al. 16 , we set A = 0.14, T e = 14 to guarantee that this signal remains subthreshold.Dξ i (t) is a Gaussian white noise with noise intensity D , satisfying Explicit Euler algorithm is utilized to solve Eq. (1) numerically, where the time step is taken as 0.001 and the initial values are set as x i (0) = 0, y i (0) = 0.

Statistical complexity measure
To well describe the collective response of neural network, we introduce the mean field of membrane potential, denoted as X(t) = 1 N N i=1 x i (t) , as the output response.Inter-spike interval (ISI) for time series X(t) is defined as the time interval between two contiguous firings of X(t) , which is recorded as mean ISI (MISI).The structural feature of MISI series is closely related to the generation of SR behavior.Here, we calculate SCM by extracting MISI time series {T s : s = 1, . . ., L} .To estimate the probability distribution from the time series of MISI, we utilize the Bandt-Pompe (BP) methodology 41 .For a given time series {T s : s = 1, . . ., L} and an embedding (1) where π = r 0 , r 1 , . . ., r (d−1) denotes the permutation of (0, 1, . . ., d − 1) .Therefore, the probability distribution P = p(π) of the corresponding ordered permutation for {T s : s = 1, . . ., L} can be written as with # being the number of corresponding permutations.
In addition, Q J [P, P e ] is the disequilibrium, quantifying the distance between the current state and the equi- librium state, which is defined as Here Q 0 is a normalization constant satisfied and J [P, P e ] represents Jensen-Shannon divergence, determined by

Quantifying noise-induced SMR behaviors by SCM
Here, we keep the model parameters ε = 0.01 , a = 1.1 , A = 0.14 , T e = 14 constant, while setting the network parameters as follows: N = 100 , k = 30 , p = 0.15 , and the coupling strength g = 0.01 .Then, we will present a discussion on the noise-induced SMR via SCM in model (1) without time delay.We will begin by extracting the MISI time series and estimating their probability distributions using the BP algorithm.Considering the requirements of the BP algorithm, it is necessary to satisfy two conditions: 3 ≤ d ≤ 7 for the embedding dimension and L >> d! for the length of the MISI time series.To ensure a more accurate evaluation of the probability distributions on the MISI time series, we set the embedding dimension to d = 3 and the length of the MISI time series to L = 60, 000 .To acquire 60,000 data points, we need to iterate through numerous output response sequences.With these settings, SCM and NSE can be obtained by Eqs. ( 4) and (5).The results are shown in Fig. 1a, b.Through careful observation, one can read that there are distinct three different noise intensities D = 0.04 , D = 0.08 and D = 0.14 , around which SCM exhibits prominent maxima, and NSE displays evident minima correspondingly.
In addition, there are extra local maximum and minimum on the curves of SCM and NSE at D = 0.001 .The low value of SCM at D = 0.001 can be attributed the distribution of the MISIs sequence approaching a uniform distribution.To highlight the local maximum of SCM and minimum of NSE at D = 0.001 , partial enlarged draw- ings of SCM and NSE with respect to noise intensity are respectively plotted in Fig. 1a, b when noise intensity varies in the range of 0.001 to 0.02.Based on these observations, we conclude that SCM achieves four maxima at D = 0.001 , D = 0.04 , D = 0.08 as well as D = 0.14 , exhibiting the striking feature of quadruple-SRs.
To verify, we plot the dependence of SNR with respect to noise intensity in Fig. 1c and present the spatiotemporal patterns of neuronal firing activities for D = 0.001 , D = 0.04 , D = 0.08 and D = 0.14 in Fig. 1d.In this analysis, we adopt the common definition of SNR, represented as SNR = P S P N , where P S is the power of the signal and P N is the power of the noise.To obtain the results for SNR, we conducted 30 numerical simulations.From Fig. 1c, it is evident that the SNR curve exhibits three maxima at nearly the same three optimal noise (2) (5) intensities as the SCM curve.However, SNR does not reach its maximum at D = 0.14 , which is different from SCM.Furthermore, Fig. 1d reveals that firing activities of neurons not only display high periodicity at D = 0.001 , D = 0.04 , D = 0.08 , but also show strong regularity at D = 0.14 , with periods approximately equal to that of the subthreshold signal.Therefore, the subthreshold signal can be detected efficiently, and both SCM and NSE attain their respective maxima and minima at these four noise intensities.We can conclude that, compared with the traditional SNR indicators, SCM and NSE indicators are more easily to detect weak periodic signals when noise intensity is strong.Thus, SCM may be a better candidate to quantify SMR behaviors and further detect the weak periodic signals in coupled excitable FHN models.Remarkably, the firing patterns of neurons exhibit noticeable differences under these four noise intensities.When D = 0.001 , neurons generate one spike in a period, while at D = 0.04 , D = 0.08 and D = 0.14 , neurons separately generate two, three and four spikes in a period.To denote these distinct firing activities, we refer to them as period-1, period-2, period-3, and period-4 firing activities, respectively.
Based on the analysis results from Li et al. 16 , we can conclude that, on one hand, the occurrence of noiseinduced quadruple-SRs can be attributed to the interaction between the internal time scale induced by the noise T 0 and the period of the subthreshold signal T e in the neuronal network.T 0 denotes the firing period of neuronal activities in the absence of external signals, corresponding to the emergence of CR in the neural network.On other hand, if the period of subthreshold signal T e falls within the range of (m × T 0 , (m + 1) × T 0 ) , noise has the capability to induce SR for m times with m being a positive integer.Following the method described by Li et al. 16 , we have determined that T 0 is approximately 3.45 for our calculations.Therefore, with T e set at 14, noise can induce four occurrences of SR in our neural network model.

Detecting delay-induced SMR behaviors by SCM
In this subsection, when keeping the model parameters ε = 0.01 , a = 1.1 , A = 0.14,T e = 14 , g = 0.01 fixed, and the network parameters N = 100 , k = 30 , p = 0.15 unchanged, we mainly use SCM to investigate SMR induced by time delay with the help of noise, and further to gain more insight into the effect of time delay on detecting the subthreshold signal in small-world neural networks.Here the embedding dimension and length of MISI time series for different time delays are always fixed as d = 3 and L = 60000 .Besides, we propose a hypothesis that if neurons do not generate firing spike under certain values of time delay, SCM and NSE are treated as 0 and 1, respectively.Based on the results of noise-induced SMR in the absence of time delay, we will pay attention to study effects of time delay on detecting the subthreshold signal under these noise intensities of D = 0.02 , 0.04, 0.065, 0.08 and 0.14, respectively.Dependence of SCM and NSE on time delay with D = 0.02 , 0.04, 0.065, 0.08 and 0.14 is shown in Fig. 2. From Fig. 2a, it is evident that for D = 0.02 , the SCM curve exhibits a series of local maxima at approximately τ =2 , 12, 26, 40, 54, 68, 82, and 96, while also achieving large values at τ =14 and 28.Conversely, the NSE curve displays an opposite trend.With the results presented in Fig. 3, we observe that neurons exhibit regular firing activities at τ =12 , 26, and 40, where their firing periods align appropriately with the values of time delay.This observation suggests a strong coherence between time delay and neuronal firing dynamics at these specific time delay values.However, at τ =14 and 28, neuronal firing activities demonstrate periodicity with their periods closely resembling  www.nature.com/scientificreports/ the period of the subthreshold signal.This finding establishes a significant correlation between the subthreshold signal and neuronal firing dynamics at these two values of time delay, leading to the occurrence of SMR.Notably, we can see that SCM does not achieve optimal values at τ = 56, 70, 84 and 98.Additionally, as evidenced by the results in Fig. 3, the neuronal firing activities at τ = 70 and 98 are non-periodic, further confirming that the subthreshold signal cannot be detected successfully and SMR does not occur in these cases.Therefore, it is summarized that under a low noise level, time delay can induce the emergence of multiple-CRs when the time delay is small or fixed around values of nT e − 2 , with n being a positive integer.Additionally, time delay can give rise to SMR occurring at nT e , where n is a small positive integer under low noise conditions.In particular, it is worth mentioning that at τ = 26 and 40, the values of SCM are slightly larger than at τ = 28 , suggesting that in the case of D = 0.02 , the strength of CR is quite greater than that of SR.How can we distinguish between SR and CR by observing the curve of SCM with respect to time delay?We can identify this based on the correlation between time delay and the period of subthreshold signal.If SCM reaches a relatively large value at an integer multiple of the period of subthreshold signal, it is considered as exhibiting SR induced by time delay in the system.On the other hand, if a local maximum is achieved at time delay deviating from an integer multiple of the period of subthreshold signal, it is believed to experience delay-induced CR in the system.Surprisingly, the investigation reveals that with an increase of noise intensity, such as D = 0.04 , 0.065, and 0.08, as presented in Fig. 2b-d, several noticeable local maxima of SCM and corresponding pronounced local minima of NSE occur at τ =14 , 28, 42, 56, 70, 84, and 98.These observations differ from the findings depicted in Fig. 2a.These outcomes, illustrated in Fig. 2b-d, indicate that delay-induced SMR phenomena can manifest at τ = nT e , where n represents a positive integer.This implies that the detection ability of subthreshold signal is enhanced when the values of time delay equal to multiples of the period of subthreshold signal, especially under moderate levels of noise.Furthermore, it is worth noting that under these noise levels, other prominent maxima may also appear on the SCM curves when time delay is small.Taking D = 0.04 and 0.065 as examples, it can be observed that optimal values of time delay, namely 8 and 2, respectively, lead to local maxima on the SCM curves.This phenomenon can be attributed to the strong correlation between time delay and neuronal firing dynamics.Interestingly, when noise intensity becomes very high, as shown in Fig. 2e for D = 0.14 , both SCM and NSE trajectories become chaotic with increasing time delay.Combining the findings presented in Fig. 2b-e, it can be concluded that under high levels of noise, the delay-induced SMR phenomena are decreased, resulting in a weakened ability to detect the subthreshold signal.Moreover, for these four noise intensities, between the optimal values of time delay, multiple small peaks can be observed on the SCM curves, suggesting the existence of multiple resonance-like phenomena between delay-induced SMR.Importantly, we have deduced that in the absence of time delay, the subthreshold signal cannot be efficiently detected at D = 0.02 and 0.065 due to the presence of local minima on SCM curves, as depicted in Fig. 1a.However, when time delay is taken into account in the information transmission between neurons, the detection capability of subthreshold signal is significantly enhanced under these two noise intensities, as demonstrated in Fig. 2c, d.Thus, we can conclude that time delay can effectively assist noise in detecting the subthreshold signal.
To further investigate the role of time delay in signal detection, we present the spatiotemporal patterns of neuronal firing activities for different time delays under D = 0.04 , as depicted in Fig. 4 positive integer (e.g., τ = 14 , 28, 84, and 98), the firing dynamics of neurons resemble those without time delay.The periods of firing activities during these time delays are nearly identical to the period of subthreshold signal.This implies that, under D = 0.04 , the subthreshold signal can be efficiently detected and the delay-induced SMR phenomena occur at τ = nT e , with n being a positive integer.However, for τ = 8 , 20, and 46, the periods of firing activities deviate from the period of the subthreshold signal, leading to a failure in detecting the subthreshold signal.Particularly, the period of firing activities at τ = 8 is almost equal to the value of time delay, indicating a strong coherence between time delay and neuronal firing dynamics.In conclusion, it can be inferred that when time delay is equal to multiples of the period of the subthreshold signal, time delay can enhance the detection ability of subthreshold signal under moderate noise levels, resulting in the delay-induced SMR phenomena.

Mechanism of delay-induced SMR
In this subsection, we will explain the phenomena of time delay-induced SMR.The results shown in Fig. 2 demonstrate that, under appropriate noise intensities, the SMR behaviors induced by time delay can occur at the integer multiples of the period of subthreshold signal.We believe that these phenomena can be attributed to the coupling effect among the time delay, the subthreshold signal and the firing dynamics of neurons.In detail, when time delay aligns with integer multiples of the subthreshold signal's period, the coupling among these factors leads to mutual activation among neurons, enhancing signal propagation efficiency and guiding neuronal firing dynamics to the optimal response state as delay-free case.As a result, neurons will exhibit synchronized firing patterns that align with the subthreshold signal's period.This synchronization greatly enhances the sensitivity of the neural network to the subthreshold signal, ultimately resulting in the occurrence of delay-induced SMR.
For example, under the noise intensity of D = 0.04 , the findings depicted in Fig. 4 reveal that neurons display synchronized and periodic firing activities when the time delay matches integer multiples of the subthreshold signal period.Furthermore, as shown in Fig. 1d, these periodic firing activities are consistent with the firing activities of neurons without time delay.Despite generating two firing activities within one period, the period of neuronal firing patterns equals the period of subthreshold signal, indicating a strong correlation between neuronal firing dynamics and the subthreshold signal during these instances.Consequently, the neural network achieves maximum perceptual sensitivity to the subthreshold signal, resulting in the emergence of multiple resonant responses at integer multiples of the subthreshold signal period.However, when the values of time delay take 8, 20, and 46, the firing patterns of neurons, although relatively regular, do not match the period of subthreshold signal, indicating no correlation between neuronal firing dynamics and the subthreshold signal.This phenomenon is likely due to the coupling effect between time delay and neuronal firing dynamics, which alters the firing pattern of neurons and gives rise to a coherent dynamic of neuronal firing.Therefore, on the SCM curves that vary with the time delay, a series of prominent peaks occur when the values of time delay equal integer multiples of the subthreshold signal's period, while smaller peaks appear at some specific values of time delay.

Effects of network parameters on delay-induced SMR
Next, we will examine the effects of various network parameters on the occurrence of time delay-induced SMR.Specifically, we focus on the rewiring probability p , degree k , and network size N , while keeping other model parameters ε = 0.01 , a = 1.1 , A = 0.14 , T e = 14 and g = 0.01 fixed.The findings are presented in Fig. 5. Analysis of Fig. 5a, b reveals that when k is fixed at 30 and N is set at 100, SCM still exhibits several maxima at τ = nT e for different values of p , whereas NSE displays an inverse pattern.It is noteworthy that the magnitudes of each peak on the SCM curve and each valley on the NSE curve do not vary significantly.This suggests that the occurrence of delay-induced SMR is robust to the rewiring probability p , and the detection ability of subthreshold signal is minimally affected by variations in the rewiring probability.
Similarly, when p is held constant at 0.15 and N remains at 100, Fig. 5c, d demonstrate that a series of maxima of SCM always occur at τ = nT e , with the increase of k , while the opposite pattern is presented on the NSE curve.This indicates the robustness of delay-induced SMR against the degree k .However, it is evident that when τ < 50 , the height of each peak on the SCM curve first increases and then decreases with the increase of k .This implies the existence of an optimal value of k that facilitates the enhancement of delay-induced SMR under moderate time delays.But, when τ > 50 , the height of every peak of SCM decreases as k increases, suggesting that higher degrees of small-world networks may weaken the strength of delay-induced SMR for large time delays.Consequently, by adjusting the degrees of small-world networks, we can effectively control the phenomena of delay-induced SMR and maximize the detectability of subthreshold signal.
Moreover, when p remains at 0.15 and k stays constant at 30, SCM and NSE display a series of maxima and minima at τ = nT e for various values of N , as illustrated in Fig. 5e, f.This highlights the robustness of delay- induced SMR in relation to the network size N .Upon careful examination, it is evident that larger network sizes result in higher peak heights for SCM and lower valley depths for NSE.This suggests that larger network sizes effectively strengthen delay-induced SMR.Therefore, to promote the detection efficiency of subthreshold signal, it is advantageous to appropriately expand the network sizes of small-world networks.

Effect of coupling strength on delay-induced SMR
When keeping the model parameters ε = 0.01 , a = 1.1 , A = 0.14 , T e = 14 constant, and the network parameters N = 100 , k = 30 , p = 0.15 unchanged, we explore the impact of coupling strength g on delay-induced SMR, as displayed in Fig. 6.From the figure, it can be observed that for τ < 50 , the peaks of SCM initially increase and then decrease with the rising g , while NSE shows the opposite trend.This suggests the presence of an optimal coupling strength that maximizes the intensity of delay-induced SMR within suitable delay values.However, for τ > 50 , the peaks of SCM decrease while the valleys of NSE increase as g increases.This indicates that weaker coupling strengths are more advantageous for enhancing the effect of delay-induced SMR when dealing with larger delay values.Hence, by adjusting the coupling strength, we can optimize the occurrence of delay-induced SMR and maximize the detection capability of subthreshold signal.

Conclusion
The current research is devoted to detecting and quantifying delay-induced SMR by virtue of SCM in smallworld neural networks, which are locally modeled by FHN neurons.MISI time series have been firstly extracted and BP algorithm has been employed to evaluate their probability distributions, based on which SCM and NSE are defined and calculated numerically.Without delay, there are four maxima of SCM and corresponding four minima of NSE at four optimal noise levels, indicating the generation of SMR phenomena induced by noise.SNR has been also carried out to characterize the noise-induced SMR.It has been found that SCM may be a better quantifier to detect SMR in neural networks than SNR.However, when time delay is introduced into the coupling process, SR can be either enhanced or decreased by time delay in small-world neural networks.At a low noise level, an intriguing phenomenon occurs: an increasing time delay leads to the emergence of regular www.nature.com/scientificreports/spatial-temporal orders in neuronal firing dynamics.Furthermore, it is observed that the periods of regular firing activities in neurons align closely with the optimal values of time delay, which indicates the strong coherence between time delay and neuronal firing dynamics.Indeed, the firing periods of neurons will undergo variations in response to changes in time delay.In the presence of moderate noise intensity, SMR phenomena can occur at integer multiples of the subthreshold signal's period, which are characterized by some prominent local maxima of SCM.Under high noise intensity conditions, SCM exhibits irregularity with varying time delay, leading to the weakness of delay-induced SMR.Finally, we have revealed the robustness of network parameters and coupling strength in the occurrence of delay-induced SMR, and demonstrated that the degrees and sizes of small-world networks, as well as the coupling strength can regulate the strength of delay-induced SMR.
It is well known that weak signals are ubiquitous in brain and usually encode the neural information.We thus expect that SCM, as an effective tool, can be used in weak signals detection in neural systems.There remains one lingering problem on how to evaluate the optimal values for subthreshold signals that are masked by noise.Perhaps we can refer to the study by Liang et al. 42 , where they controlled parameter diversity and coupling strength to jointly modulate the waveform or period of collective activity in the system, resulting in resonance for optimal suprathreshold signals.This subject holds significant importance for applications in the field of signal processing.

Figure 1 .
Figure 1.Dependence of SCM, NSE and SNR with respect to noise intensity is shown in (a-c), respectively; (d) space-time plots with noise intensity D = 0.001 , D = 0.04 , D = 0.08 and D = 0.14.

Figure 3 .
Figure 3. Space-time plots with different values of time delay under D = 0.02.

Figure 4 .
Figure 4. Space-time plots with different values of time delay under D = 0.04.

Figure 6 .
Figure 6.Dependence of SCM and NSE on different coupling strengths with respect to time delay.